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III. THE PARAXIAL WAVE EQUATION -- PROPAGATION OF GAUSSIAN 
BEAMS IN UNIFORM MEDIA 

Derivation of Paraxial Wave Equation: 

In point-to-point communication, we may think of the electromagnetic field as propagating 
in a kind of "searchlight" mode — i.e. a beam of finite width that propagates in some 
particular direction. In analyzing this mode of wave propagation, we make use of an 
important solution to the so call paraxial approximation of the electromagnetic wave 
equation (or, more precisely, the paraxial approximation of the Helmholz equation). 

To that end, we first derive the paraxial approximation and then examine the free-space 
Gaussian Beam solution(s). We start with the homogeneous Helmholz equation for 
the vector potential in the form — see Equation [ I- 13a ] 



We are looking for a wave propagating in, say, the z-direction, so we write a particular 
component of the potential in the form 



The function ^(r^co) represents a spatial modulation or "masking" of a plane wave 
propagating in the z-direction. The z-direction is obviously special and it is useful to 
appropriately parse the differential operators. For the grad operator we may write 



V 2 A(r,co)+(0 2 (i 8((o) A(r,Q))=V 2 A(f,(d) + k 2 A(r,0)) =0 



[III-l] 



A a (r,co) = v P(r,co)exp(-z'fcz) 



[III-2] 



grad { 



}=V{ 



}=v t { 




} [HI-3] 



where, for example, 
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v,{ 



}-if{ 

a x 



} + yf { 

dy 



so that 



[III-4] 



V A a (f ,co ) = ^ V, ¥ (r,co ) + z A- ¥ (r,co )-ikz r,co ) j exp(-i* z) [ HI-5 ] 



For the Laplacian operator we may write 



A a (r,co) =V 2 ¥(r» exp(-/£ z ) + A |^ A <i>(r,co) - /Jk ¥(r,co) j exp(-iikz)j 



[III-6] 



where, for example, 

v t 2 { 

Therefore, 



i 32 ; 



V 2 A a (r,co)=j^ V 2 ¥(r» + ^ ¥(r,a>)-2ifc |- ¥(r,co)-fc 2 ¥(r,co)j exp(-zfcz) 



[III-7] 



[III-8] 



and the parsed Helmholz equation ( without approximation ) becomes 

V t 2 Y( r, to) - 2 /it — Y (r,co ) + — - Y (r,co ) = 

OZ oz 



[III-9] 



The paraxial approximation is precisely defined by the condition 

d 2 



_d_ 

dz 



2ik — ¥(?,©) » — - ¥(?,©) 

oz 



[III- 10] 
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which means that the longitudinal variation in the modulation function, ^(r^co ) , changes 
very little in the wavelength associated with beam ~ i.e. 2n/k . In this approximation, we 
neglect the third term and obtain the equation 



V t 2 ¥(r,co)- 2ik — ¥(?,©) = 

oz 



[III-ll] 



which is called the paraxial approximation of the wave equation. 15 

Solutions of the Paraxial Wave Equation 
The Gaussian beam 

To inform or motivate our next step, we consider the paraxial approximation of a 
known solution of the Helmholz equation ~ i. e. a spherical wave 



exp 



(-ikr) _ exp(-i* Jx 2 + y 2 +z 2 ) 



exp 



Jx 2 + y 2 +z 2 




exp(-ikz )exp 



-ik [x 2 + y 2 ) 



2z 



[111-12] 



Reflecting on the "quadratic" form of this approximate expression, it is reasonable to 
look for an axially symmetric solution of the paraxial wave equation in the following 
form - i.e. a Gaussian beam: 



Y G (p , z, 0) ) = A G exp[ -i P(z)] exp 



ik p 2 
2q(T) 



[111-13 ] 



Obvious the paraxial equation has the same mathematical form as the Schrodinger equation and, thus, all that is 
know about solutions of that equation may be directly applied to understand issues in light propagation (or visa 
versa). 
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where p 2 =x 2 + y 2 . 

We test our conjecture by substituting the Gaussian beam function — i.e. Equation 
[ III- 13 ] ~ into the paraxial wave equation — i.e. Equation [ III- 1 1 ] ~ to wit 

exp H p (Z )] v? {-fSF 2ik I. l exph ' p w] exp r3> 1 =0 • [ m - 14 1 



Executing the indicated operations, we obtain 



r ikp 2 i\ nk k 2 p 2 „. r . a ikp 2 b Ti 



expHP( z) ]exp r — jj-— _ 



-2/Jt 



- / ^ P ^)+ ,;' V ^ 1 2 ^<?(z) ^ = [III-15a] 
3z 2[^(z)] 3z Jj 



or simplifying 



„, a „, , nk k 2 p 2 f a , A n 



[III-l5b] 



Hence, for an arbitrary p this equation is separable into two parts ~ viz. 



— q(z)=l [III-15C] 
a z 



and 



2fc— P(z) + — =0 
oz q{z) 



^ P ( ) _ _ Z 
3 z q{z) 



[III-l5d] 



which are satisfied by the simple solutions 

q(z) = z+q 



[III-l6a] 
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and 

— P( z ) = ^- = -i — ln[z + q ] ^ P(z) = -iln[z+? ] . [III-16b] 

az Z+q az 

On comparison with the paraxial approximation of a spherical wave — i.e Equation 
[ III- 12 ] — we may write q{z) in terms of a radius of curvature R (z) and a width 
w(z) -- viz. 

1 1 1 -i 2 

+ 2~TT • [111-17 ] 



q(z) z + q R(z) kw 2 (z) 

To standardize the constants of integration we assume a plane wavefront at an 
arbitrary reference point z= — i.e. we take P(0) = °°. Thus, 

! =0 [III-18a] 



*(0) 

-i2 1 ikw 2 (0) iTiw 2 (0) rTTT1Q , 1 

and — ^— =— => q = = =i L F [III-18b] 

k w yO) q 2 A 

where L^ = k w 2 (0)/2 = 7t w (0)/a is the critical Gaussian beam scaling parameter 
which is called variously the Fresnel length, the diffraction length, or the 
confocal parameter. In terms of this parameter, Equation [ III- 17 ] may be written 

— ~ + -=^=^ = ^ [ni-19] 



q(z) R(z) kw 2 (z) z + iLp z +4 

Equating real and imaginary parts, we obtain 

1 _ z a ^ -i 2 -z'L F 

rJz) ~ z 2 +L 2 F ' UK kw 2 (z) ~ z 2 +L 2 F 
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R(z) = z[l + Ll/z 2 } 

w 2 (z) = w 2 (0)[l + z 2 /4] 

where Lp = % w 2 (0)/X 



[ 111-20 ] 



Now since Equation [ III- 16b ] may be written 

P(z) = - i In[z + q ] = -i ln[ z + i Lp] = -i {in [z 2 + L\] + i tan" 1 [ip/z]}/ [ III- 16b' ] 



we may write 



exp[-i P(z)] = 



exp [-z tan 1 [ /z ]] exp [- i tan 1 [L F / z] ] 



to obtain the usual, officially approved form of the Gaussian Beam 
%(p,z,a>) = i4 exp[-iP(z)]expj^-^-^j 



2 1 



. w(0) r ., _i rr . fl r ikp 2 1 f p 
= 4 — exp[-ua„ [v^—j^- j 



[ 111-21 ] 
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The following kind of picture is sometimes found to be helpful in understanding the 
propagation of a Gaussian Beam (the bold curve depicts the spatial variation of the 
beam width and the light curve the beam curvature at particular points in space): 



4.T 



w(z) 

W(0) 



-3. 



-2. 



-1. 



1 . 



2. 




w(z) = A 

Z It-Vv'(O) 



Higher order Hermite-Gaussian beams 

In order to study the propagation of higher order beams, we substitute the following 
trial solution: 



*¥ H _ G (p,z,0)) =F(x,y, z)V G (p,z,(D) 

= f g(j ) exp[-/ * (z)] %(p, z, CO) 



[ 111-22 ] 



into the paraxial wave equation - viz. Equation [ III- 11 ] - and obtain 
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F(x, y, z) V t 2 Y G (p, z, co) + 2[v t F(x, y, z) ■ V t Y G ( p, z, ©)] + V F G ( p, z, co) V 2 F(* y, z) 



-2ik^ G (p, z,(0)^- F(x,y, z)-2ik F(x,y,z)^- %(p,z,co)=0 
az a z 



[ 111-23 ] 



Since the sum of the first and fifth terms already satisfies the paraxial wave equation, 
Equation [ 111-23 ] reduces to 



[ dw w~\ f g" \ dw w~\ g' 



d<P 



— +2ik\ — j^f 7 " 1 \-2ik\— \y 2kw = 



L dz q 



g 



i dz q J 



dz 



[ 111-24 ] 



From Equations [ III- 19 ] and [ 111-20 ] we see that 



dw w _w 
dz q ~ R 



w -i 2 
R kw 



1 2 
k w 



so that the reduced equation — i.e. Equation [ 111-24 ] ~ becomes 

^_4^ + JL_ 4q JL_ 2 ^^ = o 



g 



dz 



[ 111-25 ] 



where £ = x/w and q = y/w . 



A Hermite polynomial of order n 16 has the following differential equation: 



^H„(x)-2x^-H„(x)+2nH„(x)=0. 



dx 



dx 



[ 111-26 ] 



d" 

16 The Hermite polynomials have the generator H n (x ) = (-1)" expfx 2 ) exp(-X 2 ) . 

dx 
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With the simple change in variables x = J2t l = -Jl x/w and a = -Jl q = -J2 y/w 
Equation [ 111-25 ] may be written 



l[d 2 f _ df~\ 1 [d 2 g _ dg~\ 



- J T _ 2X 1+ - | T _ZO I 

f[_dx JxJ gldo do] 



2a — i - 2 fcw 



J3> 



= 



[ 111-27 ] 



Thus, it is apparent that we can write the functions f^~~^ an d as Hermite 

polynomials ~ viz. 



f(-) = H„(x) = H„(y2^) and g( ± ) = H» = H J Jl ± 



2 



2 d<S> 



if we require that 2kw = -2(n+ m). Hence 

dz 

d<t> _ (n+m) _ (n+m)Lp 
dz kw 2 2[Lp + z 2 ] 

or $(z) = -(n + m) atan(z/L F ). 

Finally we may write a general solution for the paraxial equation as 



XI (P, z, co) = AZ ^ Hj 72 * ) Hj 72 1 



x exp[z [n + m + 1] tan ^z/^f)] ex P 



ik p 2 




P 2 1 


L 2tf( z )_ 


exp 


. w 2 (z)_ 



[ 111-28 ] 
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A GALLERY OF HERMITE-GAUSSIAN FIELD DISTRIBUTIONS 
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Gaussian Beam Transformation Matrices 

What we have shown above is that a given Hermite-Gaussian beam is essentially 
completely specified or defined by the complex function q{z) . In propagating through an 
optical system, the beams are transformed by various optical components. The amazing 
fact is that the transformation produced by a given component follows a 
simple ABCD transformation law - viz. 



[ 111-29 ] 




where A, B, C, D are the matrix elements found in our analysis of 
geometric optics!! 

To "prove" this, we argue by example. For example, the transformation through a 
uniform dielectric region is given by 

q 2 = q { + L 

so that {A = 1 ; B = L ; C = 0; D= 1} and the transformation through a thin lens is given 
by 

_L ^ _L_I 

<?2 <7l f 

so that | A = 1 ; 5=0; C= ~J' D=1 }- 



Further "justification" of this transformation law may be found in terms of the so called 
" p/p ' argument"- viz 
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■w-tgirt o f sphtrical 
wavt associated 
with ajvtn input ray 




■w-ig/rtei sphtrical 
wavt ass ociattd with 
givtn output ray 



From geometric optics and, in particular, Equation [ II- 1 1 ], we may write 



Pout 
Pout 



\ 



A — + B 

Pin 



V 1 



C^+D 



or 



A^ t =(AAz in +B)(CAz m +D)- 1 



which is identical to transformation equation [ III- 29 ] if we interpret q(z) as the wave 
optics generalization of A z = p/ p ' . 
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